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ABSTRACT 

We apply the Jeans equation to estimate masses for eight of the brightest dwarf spheroidal (dSph) 
galaxies. For Fornax, the dSph with the largest kinematic data set, we obtain a model-independent 
constraint on the maximum-circular velocity, Vmax = ■ Although we obtain only lower- 

limits of Vmax ^ 10 km for the remaining dSphs, we find that in all cases the enclosed mass at the 
projected half-light radius is well constrained and robust to a wide range of halo models and velocity 
anisotropics. We derive a simple analytic formula that estimates M{rhaif) accurately with respect to 
results from the full Jeans analysis. Applying this formula to the entire population of Local Group 
dSphs with published kinematic data, we demonstrate a correlation such that M{rfiaif) oc "^^haif^^j O'' 
in terms of the mean density interior to the half-light radius, (p) oc "r^aij^*^ '^- This relation is driven 
by the fact that the dSph data exhibit a correlation between global velocity dispersion and half-light 
radius. We argue that tidal forces are unlikely to have introduced this relation, but tides may have 
increased the scatter and/or altered the slope. While the data are well described by mass profiles 
ranging over a factor of < 2 in normalization {Vmax ~ 10 — 20 km s""'^), we consider the hypothesis 
that all dSphs are embedded within a "universal" dark matter halo. We show that in addition to the 
power law M cx r^'^, viable candidates include a cuspy "NFW" halo with Vmax ^ 15 km s~^ and 
scale radius tq ~ 800 pc, as well as a cored halo with Vmax ~ 13 km s~^ and tq ~ 150 pc. Finally, 
assuming that their measured velocity dispersions accurately reflect their masses, the smallest dSphs 
now allow us to resolve dSph densities at radii as small as a few tens of pc. At these small scales we 
find mean densities as large as (p) < 5M0pc~'^ (< 200GeV cm~'^). 

Subject headings: galaxies: dwarf — galaxies: kinematics and dynamics — (galaxies:) Local Group 
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1. INTRODUCTION 

Dwarf spheroidal (dSph) galaxies are the smallest and 
least luminous galaxies in the Universe, and so provide 
unique diagnostics of ga laxy formation at small scales. 
Early observations (e.g.. IAaronsonlfl983f ) of stellar kine- 
matics in the brightest {Ly ~ 10^~'' Lv,q) of the Milky 
Way's dSph satellites yielded the surprising result that 
all well-studied dSphs have central velocity dispersions of 
~ 10 km s~^, larger than expected for self-gravitating, 
equilibrium stellar populations with scale radii of ~ 100 
pc. Following work showing that such dispersions are 
unlikely to hav e been inflated signifi c antly by uniden- 
tified binaries (jOlszewski et al.l Il996t iHargreaves et al 
1996D . stellar- atmospheric turb ulence (e. g., Prvo r et al 



1988!) or tidal disruption (Piatek fc Prvoiiil9 95: Oh et al 



1995). it has become widely accepted that dSphs are 
dominated by dar k matter, wi th mass-to-light ratios of 
M/Lv ~ 10^"^ (jMateol [l998l and references therein). 
Furthermore, that dSphs have similar velocity disper- 
sions despite varying over several orders of magnitude 
in luminosity suggests that the dark matter halos of 
dSphs are more similar than the differences in their stel- 
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lar c omponents rnight o therwise suggest. Along these 
lines. [Mateo et ahl (|1993D noted that all dSphs measured 
at that time had dynamical masses of ^ IO^Mq, im- 
plying an anti-correlation between luminosity and M/L 
(see, e.g.. Figure 9 of,Mateo 1998). The recent disc overies 
of ultra-faint Milky Way s at ellites (e.g., [W illman et ahl 
I2005al: IZucker et all I2006bl lal: IBelokurov et al.., 20071 ) ex- 
tend the range of dSph structural parameters by an or- 
der of magnitude in (luminous) scale radius and by three 
orders of magnitude in luminosity. These objects pro- 
vide new opportunities to identify scaling relations that 
may point to unifying principles tha t gove rn dSph for- 
mation and evolution. IStrigari et "all (|2008l . "SOS" here- 
after) present one such (non-) rel ation, sharpening the 
point made bv iMateo et al.l (|1993D by arguing that, de- 
spite spanning five orders of magnitude in luminosity, all 
the studied dSphs have mass ~ 1O''-M0 enclosed within 
their central 300 pc. S08 contend that a common value 
of M300 ~ IO^Mq may represent the minimum mass of 
a dark matter halo, or it may represent the minimum 
mass required for galaxy formation (e.g.. iLi et al.|[2009l : 
Maccio et al. 2009). 

Here we present an analysis that differs from that 
of S08 in one key respect. Rather than M300, we 
consider the mass enclosed at the half-light radius, 
M{rhaif), which is well constrained by the available kine- 
matic data. Our reasons for this choice all stem from 
the fact that, among the Local Group's known dSphs, 
fhai f varies over two ord ers of magnitude, from a few 
tens of pc to ~ 1 kpc (jlrwin fc Hatzidimitrioul 119951 : 
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iMcConnachie fc Irwinl 120061 : iMartin et al.ll2008f) . First, 
regarding the faintest dSphs for which all detected mem- 
bers lie within projected radii R < 100 pc, M300 has 
literal meaning only if one extrapolates the properties of 
the inferred dark matter halo far beyond regions sampled 
by observational data. Since there is no evidence that the 
dark matter halos of the smallest systems extend to radii 
of order 300 pc, for these objects M300 is merely a param- 
eter setting the normalization of the mass profile (Sec- 
tion l3.3p . Second, a common value of -M300 among dSphs 
relays little information about their total masses — e.g., 
one can build a plausible Milky Way model using a cen- 
tral halo with V^ax ~ 200 km s"^ and M300 ~ lO^M©. 
Third, for the smallest galaxies, constraints on M{rhaif ) 
imply constraints on dark matter densities in regions as 
small as r 30 pc, providing unprecedented mass resolu- 
tion on small galactic scales. Finally, by considering the 
empirical relation we derive between M{rhaif) and rhaif, 
we find reason to make a stronger claim than SOS's con- 
cerning the common mass of dSph satellites. Specifically, 
we find that the entire population of Local Group dSphs 
can be fit reasonably well with a "universal" dark matter 
halo, about which the scatter with respect to the data is 
similar to that about a common value of M300. That is, 
to the extent that dSphs can be said to have a common 
M300, the data suggest dSphs may have common masses 
at all radii. 

2. DATA 

The analysis presented below proceeds in two stages. 
In the first (Section [3]), we apply the Jeans equation to 
estimate halo parameters of the eight dSphs for which we 
have acquired large kinematic data sets. For the Carina, 
Fornax, Sculptor and Sextans dSphs we use kinematic 
data obtained with the Michi gan/MIKE Fi b er Spe ctro- 
graph (MMFS) at Magellan. I Walker et all (|2007aD de- 
scribe the acquisition and reduction of these data, while 
iWalker et al.l ()2009bD present the entire data set. For 
Draco, Leo I, Leo II and Ursa Minor we use data gathered 
with the He c toche lle fiber spectrograph at the MMT. 
iMateo et al.l ([2008') describe the acquisition and reduc- 
tion of these data and present the data for Leo I. Pa- 
pers presenting Hectochelle data for the remaining three 
galaxies are in preparation. 

One purpose of our Jeans analysis is to justify the use 
of a simpler method of estimating enclosed masses (Sec- 
tion 13. 6p , which we can then apply to a broader sample 
of objects. In the second stage of this work (Section 
|3|), we demonstrate correlations among the bulk struc- 
tural and kinematic properties of dSphs. There we ben- 
efit by considering the entire dSph population, includ- 
ing not only the bright Milky Way satellites b ut also 
the ultra-faint satellites disc overed (e.g., Willma n et al.l 
I2005bl : iZucker et al. 2006b. ai: iBelokurov et al...200^ and 
observed spectroscopically (e.g., IMartin et al.l 120071 : 
ISimon fc Gehall2007l: IBelokurov et al.l 120091) during the 
past five years. We also include three dSph satellites of 
M31 (Andll, AndIX and AndXV) that have published 
kinematic dat a (we do not include AndXVI, for which 
(|Letarte et al.l 2009) report an unresolved velocity dis- 
persion), as well as the two remote dSphs (Tucana and 
Cetus) at the Local Group's outskirts. For all of these 
objects. Table |T] lists (projected) half-light radii, velocity 
dispersions and luminosities that we adopt from either 



the literature or our own data. 

In the interest of uniformity, we have chosen where 
possible to adopt measurements from studies of multiple 
objects observed and analyzed using the same method- 
ology. We caution that in some cases this preference 
may introduce systematic bias. For example, for most 
of the low-luminosity dSphs we adopt t he structural 
parameters measured from SDSS data by IMartin et al.l 
(2008). In our experience, deeper follow-up imaging for 
the faintest objects often indicates smaller half-light radii 
than are apparent from SDSS data alone — compare, for 
example, the stellar maps of Segue 2 from SPSS and 
MMT/Megacam imaging in Figure 3 of IBelokurov et aLl 
( 20_ 0i^ (see also the LBT follow-up imaging of Hercules 
bv ISand et al.l l2009f) . In most cases the difference is 
slight, but we caution that the structural parameters of 
the faintest systems may be subject to revision when 
deeper data become available for more of these objects. 

2.1. Velocity Dispersion Profiles 

For the Jeans analysis in Section [3] we use empiri- 
cal, projected velocity dispersion profiles fo r eight bright 
dSphs with abundant kinematic data. IWalker et al.l 
{2007M) present velocity dispersion profiles for seven of 
these objects, but these profiles can now be updated — 
and we can now provide a profile for Ursa Minor — using 
the data we have obtained during the past two years. For 
each galaxy we compute the velocity dispersion profile 
after 1) discarding all stars for which the probability of 
dS ph membership, acco rding to the algorithm described 
bv IWalker eFall (|2009cD . is less than 0.95; 2) subtract- 
ing the mild velocity gradients t hat likely refiect the bulk 
transverse motion of the dSph (IWalker et al.ll2008l ). and 
binning the data using circular annuli containing equal 
numbers of member stars. We estimate the velocity dis- 
persion in each bin using the maximum -likelihood tech- 
nique described bv IWalker eFaLl (|2006D . 

Figure [Tj displays the velocity dispersion profiles we 
measure for Carina, Draco, Fornax, Leo I, Leo II, Sculp- 
tor, Sextans and Ursa Minor. Whi le th e profiles remain 
generally flat (c.f. Figure 2 of Walker et al. 2007b), we 
detect a gentle decline in velocity dispersion at large 
projected radii (i? > 1 kpc) in Fornax. This behav- 
ior becomes appa rent only after addit ion of the most re- 
cent MMFS data (jWalker et al.ll2009H l. which contribute 
most of the data points at large radius. We also find hints 
of gently declining profiles in Sculptor and Sextans, al- 
though these data are noisier and the apparent declines 
hinge on the outermost o ne or two da t a poin ts. 

Usin g the same data, iLokas et"all (|2008l ) and iLokad 
(|2009l ) have recently measured profiles for Carina, For- 
nax, Leo I, Sculptor and Sextans that all decline more 
significantly at large radii than the profiles we present 
in Figure [Tj The d ifferen ce arises because ILokas et al.l 
(l2008i) and iLokasI (j2009f ) adopt a rejection algorithm 
that attempts to remove stars that may have been 
tidally stripped and thus do not provide reliable trac- 
ers of the gravitati onal potential. The algorithm, de- 
scribed in detai l by Kliment owski et al. (2007) (see also 
Iden Hartog fc Katgert T996i) , iterativcly estimates the 
mass profile from the velocities and positions of accepted 
members and then rejects stars with velocities that dif- 
fer from the mean by more than ^^2GAI(r)/r. The va- 
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TABLE 1 

dSph Structural Parameters, Velocity Dispersions and Estimated Masses* 



Object Lv Thaif o-Vq Af(fhalf) (p> Ref. 



Carina 


2.4 ±1.0 


X 


10"" 


/I 1 1 00 

241 ± 26 


6 


6 ± 1.2 


6.1 ± 2.3 


X 


10 


1.0 ± 


4 


X 


1 n— 1 
10 


1,2 


Draco 


2.7 it 0.4 


X 


lO'' 


196 zb 12 


9 


1 it 1.2 


9.4 it 2.5 


X 


iU 


3.0 ± 


8 


X 


1 n— 1 
iO 


3,4 


Fornax 


1.4 it 0.4 


X 


iO 


C^a Q _L Q A 
DDO ± o4 


11 


7 it 0.9 


5.3 it 0.9 


X 


10 


4.2 ± 


7 


X 


10 


1,2 


Leo I 


3.4 it 1.1 


X 


10^ 


246 lb 19 


9 


2 it 1.4 


1.2 it 0.4 


X 


1 n7 
10 


1.9 ± 


6 


X 


1 n— 1 
10 


1,5 


Leo II 


5.9 it 1.8 


X 


10^ 


151 lb 17 


6 


6 it 0.7 


3.8 it 0.9 


X 


10^ 


2.6 ± 


6 


X 


1 n— 1 
10 


1,6 


Dculptor 


1.4 it 0.6 


X 


10^ 


260 ± 39 


9 


2 it 1.1 


1.3 it 0.4 


X 


1 n7 
10' 


1.7 ± 


5 


X 


T n— 1 
10 


1,2 


Sextans 


4.1 it 1.9 


X 


10^ 


682 it 117 


7 


9 it 1.3 


2.5 ±0.9 


X 


10^ 


1.9 ±0 


7 


X 


10^2 


1,2 


UMi 


2.0 it 0.9 


X 


10^ 


280 it 15 


9 


5 it 1.2 


1.5 ±0.4 


X 


10^ 


1.6 ±0 


4 


X 


10-1 


1,7 


Bootes 1 


3.0 it 0.6 


X 


10 


/1 _i_ 1 

z4z it 2i 


6 


5 it 2.0 


5.9 ± 3.7 


X 


1 n6 
10 


1.0 ± 


6 


X 


1 n— 1 
10 


3,8 


riootes Z 


1.0 it 0.8 


X 


10^ 


51 it 17 


10 


5 it 7.4 


3.3 ± 3.3 


X 


10^ 




5 


9 


± 5.9 


3,9 


UVen 1 


2.3 it 0.3 


X 


10 


t^tiA _1_ Oti 

oo4 ± 3o 


7 


6 it 0.4 


1.9 ± 0.2 


X 


10 


2.5 ± 


3 


X 


10—2 


3,10 


f^\fan TT 

ven 11 


( .y It 0.0 


X 




/4 It iZ 




D It i.U 


y . i It 4.Z 


X 


iU 


0.0 It Z 





X 


in-i 


o,lU 


C^oma 


3.7 it 1.7 




10^ 


77 it 10 


4 


6 it 0.8 


9.4 ± 3.5 




10^ 


4.9 ± 1 


8 




10-1 


3 10 


Hercules 


3.6 it 1.1 


X 


10^ 


330 it 63 


3 


7 it 0.9 


2.6 ± 1.4 


X 


10^ 


1.7±0 


9 


X 


10-2 


3^1 


Leo IV 


8.7 it 4.6 


X 


10^ 


116 it 30 


3 


3it 1.7 


7.3 ±7.3 


X 


10^ 


1.1 ±1 


1 


X 


10-1 


3,10 


Leo V 


4.5 it 2.6 


X 


10^ 


42 it 5 


2 


4 it 1.9 


1.4 ± 1.4 


X 


10^ 


4.5 ±4 


5 


X 


10-1 


12,13 


Leo T 


5.9 it 1.8 


X 


10* 


178 it 39 


7 


5 it 1.6 


5.8 ±2.8 


X 


10'^ 


2.5 ± 1 


2 


X 


10-1 


3,10,14 


Segue 1 


3.3 it 2.1 


X 


102 


29 it 7 


4 


3 it 1.2 


3.1 ± 1.9 


X 


10^ 




3 





± 1.8 


3,15 


Segue 2 


8.5 it 1.7 


X 


10^ 


34 it 5 


3 


4 it 1.8 


2.3 ±2.3 


X 


10^ 




1 


3 


± 1.3 


16 


UMa I 


1.4 it 0.4 


X 


10^ 


318 it 45 


11 


9 it 3.5 


2.6 ± 1.6 


X 


10^ 


2.0 ± 1 


2 


X 


10-1 


3,8 


UMa II 


4.0 it 1.9 


X 


10^ 


140 it 25 


6 


7it 1.4 


3.6 ± 1.6 


X 


10^ 


3.2 ± 1 


4 


X 


10-1 


3,10 


Willman 1 


1.0 it 0.7 


X 


10^ 


25 it 6 


4 


3 it 1.8 


2.7 ±2.3 


X 


10^ 




4 


1 


±3.6 


3,8 


Andll 


9.3 it 2.0 


X 


10** 


1230 it 20 


9 


3 it 2.7 


6.2 ±3.6 


X 


10^ 


7.9 ±4 


5 


X 


10-3 


17,18 


AndIX 


1.8 it 0.4 


X 


10^ 


530 it 110 


6 


8 it 2.5 


1.4± 1.1 


X 


10^ 


2.3 ± 1 


7 


X 


10-2 


19 


AndXV 


7.1 it 1.4 


X 


10^ 


270 it 30 




iiite 


1.9 ± 1.9 


X 


107 


2.3 ±2 


3 


X 


10-1 


20,21 


Cetus 


2.8 it 0.9 


X 


lO'' 


590 it 20 




17 it 2 


9.9 ±2.3 


X 


10^ 


1.1±0 


2 


X 


10-1 


17,22 


Sgr*** 


1.7 it 0.3 


X 


10^ 


1550 it 50 


11 


4 it 0.7 


1.2 ±0.1 


X 


10« 


7.5 ± 1 





X 


10-3 


23,24 


Tucana 


5.6 it 1.6 


X 


10^ 


270 it 40 


15 


8 it 3.6 


4.0 ± 1.9 


X 


10^ 


4.7 ±2 


3 


X 


10-1 


25,26 



Estimated using Equation 11 
Reference s : 1) llrwin & : Hatzidimitriou' figgSj) ; 2) IWal kcr ct al.' f2009dP; 3) IMartin etldl I I2008D : 4) 
IWallcer et"all Il2007bl); 5) fMatco ct al. (2008); 6) IKocli e~ al. (2007a): 7) Walke r et al. in prepar ation; 
8) 'Martin ct al.' (2007); 9) Kocli ct al. (2009); 10) 'STinon & Gcha (2007); 11) Aden ct al." (20091): 12) 
[Bclokurov ct al. (2008); 13) Walker ct al. (2009a); 14) Irwin ct al. (2007); 15) Gcha ct al. (20 09al); 16) 
[Bclokurov ct al. (2009): 17) McConnachic & Irwin (2006); 18) Cote ct al. (1999i') : 19) iChapman et al.l 
(2005); 20) Ibata ct al. (2007); 21) Lctartc ct al. (2009); 22) Lewis ct al. (2007); 23) llbata &: IrwinI (119971) ; 
24) jMaicwski ct al., (.2003i ): 25) iSaviane ct al.. tl99a) ; 26) iFratcrnali ct al.. {200% 

Structural parameters refer to the bound central region of Sgr (see lMaiewski et al.ll2003l) . 



lidity of this method clearly depends on the accuracy 
of the mass pro file estimate , and this q uality is diffi - 
cult to assess. ILokas et all (|2008f ) and ILokasI H2009f ) 
use a mass estimator derived from the virial theorem 
IJleislcr ct al. 1985), which holds when the tracer popula- 
tion is self-gravitating or, if there is a dark matter compo- 
nent, when the mass-to- light ratio is constant with radius 
("mass follows light"). Then it is not surprising that af - 
ter i mposing this rejection algorithm ILokas et al.l (120081 ) 
and ILokasI (j2009D obtain falling velocity dispersion pro- 
files t hat are consistent wi th mass-foUows-light models. 
While iKlimentowski et al.l (|2007) succesfuUy test their 
algorithm and mass esti mator using an N-bq dy simula- 
tion of a tidally "stirred" (|Maver et al.ll2001allbf i satellite 
in which mass actually does follow light (after 99% of the 
initial dark mass is stripped), they do not demonstrate 
that their method preserves flat velocity dispersion pro- 
files where they may genuinely exist — e.g., in equilibrium 



systems with extended and intact dark matter halos^. 

Pending evidence that the t idally stirred N-body model 
of iKhmentowski et al.l ()2007f ) describes the evolution of 
the real dSphs in our sample, wc opt not to apply the re- 
jection method used by Lokas ct al. (2008); Lokas (2009). 
Thus we avoid imposing any bias toward mass-foUows- 
light models that may result from use of the virial the- 
orem mass estimator, but at the cost of leaving our 
samples vulnerable to contamination from tidally un- 
bound stars. Given the difficulty of identifying such 
stars, we choose to err on the side of including them 
in our analysis. For at least some systems — e.g., Fornax, 
which has o rbital pericenter near its present distan ce of 
- 140 kpc (jPiatek et all [20071 : iWalker et all [2001 — we 
expect a negligible contribution from unbound tidal de- 

"^ IKlimentowski et al.l (120071 ) note that their method rejects fewer 
than 1% of stars drawn from an equilibrium system with an ex- 
tended dark matter halo, but they do not describe the locations of 
the rejected stars or the impact on the measured velocity dispersion 
profile. 
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bris. Among the remaining dSphs in our sample, there 
exists kinematic evidence of tidal streaming motion in 
the outer parts of Carina (iMurioz et al.ll2006l ) and Leo I 
(|Sohn et al.ll2007t iMateo et al.ll2008f l. In fitting models 
to the velocity dispersion profiles of these two galaxies, 
as well as Draco, which shows a rising profile in its outer 
parts, we experimented with discarding the outermost 
1-2 data points in each system. Since the results were 
indistinguishable from fits that included the outermost 
points, we decided to include the outer data points in 
our final analysis. 

3. JEANS ANALYSIS 

With the goal of measuring dSph masses, we assume 
that the data sample in each dSph a single, pressure- 
supported stellar population that is in dynamic equi- 
librium and traces an underlying gravitational potential 
dominated by dark matter. Further assuming spherical 
symmetry, the mass profile, M{r), of the dark matter 
halo relates to (moments of) the stellar distribution func- 
tion via the Jeans equation: 



V dr r 



GM(r) 



(1) 



where v{r), and I3{r) = 1 — v^jvl describe the 

3-dimensional density, radial velocity dispersion, and or- 
bital anisotropy, respectively, of the stellar component. 
For the special case of consta nt anisotropy, the Jea ns 
equation has the solution (e.g.. iMamon &: Lokasl[2005f) 



,"2 ^ 



^2/3-2 



v{s)M{s)ds. 



(2) 



Projecting along the line of sight, the mass profile re- 
lates to observable profiles, the projected stellar den- 
sity, /(i^), and velocity di spersion, CTp(i?), according to 
(jBinnev fc Tremaind [20081 "BT08" hereafter) 



al(R) 



1-/5 



i?2 



x/r2 - i?2 



dr. (3) 



To estimate dSph masses via the Jeans equation we 
therefore employ the following strategy: 1) adopt a sim- 
ple analytic profile for I{R) from the literature; 2) adopt 
a parametric model for M{r); and 3) find the halo pa- 
rameters that, via equations [5] and [31 best reproduce the 
empirical velocity dispersion profiles shown in Figure [T] 

3.1. Stellar Density 

Stell ar surface densi ties of dSphs are ty pically 

fit by IPlummeri ([ToIl . \Kin^ ([l96l and/or ISersid 
(119681), profiles (e.g.. Ilrwin fc Hatzidimitrioul I199F 



IMcConnachie k Irwini[2"006l : IBelokurov et al.ll2007l ). The 
Plummer profile, I{R) = i(^rL,;)-Ml + ^V^L/]"' 
where L is the total luminosity, is the simplest as it has 
only a single shape parameter, the projected half-light 
radius^. It is also the only profile with published param- 
eters for all dSphs, since the concentration parameters 
of King and Sersic profiles arc not well-constrained by 

^ Throughout this work we define r^aif as the two-dimensional 
half-hght radius — i.e., the radius of the cyUnder that eneloses half 
of the total luminosity. 



the sparse data available for the faintest dSphs. There- 
fore, in what follows we adopt the Plummer profile to 
characterize dSph stellar densities. 

Given a model I{R) for the projected stellar density, 
one recovers the 3-dimensional density from (BT08) 



v{r) 



dl dR 



dR Vi?2 - r2 



(4) 



Thus for the Plummer profile, we have i^(r) ~ 

3i(47rrL/)-Ml + ^V^L/]-'/'- 

We note that even though dSph surface brightness data 
can be fit adequately by a variet y of density profiles, 
the ch oice of profile is not trivial. lEvans. An fc Walkeii 
(|2009l ) demonstrate that, even when the gravitational po- 
tential is dominated by dark matter, the adopted shape 
of the stellar density profile can profoundly affect the 
inferred shape of M(r) at small radii. In what follows, 
while for simplicity we present only the results obtained 
using the Plummer profile, we explicitly identify any re- 
sults that are strongly sensitive to this choice (Section 

MM- 

3.2. Halo Model 

For the dark matter halo wc follow S08 (also 
iKoch et aL, , 2007b, a) in adopting a generahze d Hernquist 
profile given by (|Hernauistlll990l : IZhaolll996D 



p(r) 



Pol — 

Vo 



1 



(5) 



where the parameter a controls the sharpness of the tran- 
sition from inner slope liuir^o dlTi{p)/d\n{r) cx —7 to 
outer slope lim^— »oo 'ilii(yo)/(iln(r) cx —3. In terms of 
these parameters, the mass profile is 



M(r) 



An I s p{s)ds — 
lo 

'3-7 3 



2F1 



47rporg f r 

3-7 
-7 3 



3-7 



ro, 

7 + a 



(6) 



a 



where 2^1 (a, b', c; z) is Gauss's hypergeometric function. 

The profiles admitted by Equation [5| include the range 
of plausible halo shapes relevant to the ongoing con- 
troversy regarding "cores" versus "cusps" in individ- 
ual dark matter halos. Profiles with 7 > are cen- 
trally cusped while those with 7 = have constant- 
density cores. For a = 7 = 1, one recovers the cuspy 
iNavarro. Frenk fc White! (fT99l [T997, "NFW" hereafter) 
profile motivated by cosmological N-body simulations. 

3.3. Normalization 

It is convenient to normalize the mass profile by M(ro), 
the enclosed mass at the scale radius of the dark mat- 
ter halo. This quantity relates to the maximum circular 
velocity, Vmax , according to 



rp [ 3-7 3—7 . 3— 7+a . -i 1 
2^1 L Q ' Q ' a 'J 



3—7 3 — 7 . 3 — 74-Q . 



where r] = rmax/^o identifies the radius corresponding to 
the maximum circular velocity, and is specified uniquely 
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Fig. 1. — Projected velocity dispersion profiles for eight bright dSphs, from Magellan/MMFS and MMT/Hcctochcllc data. Over-plotted are 
profiles calculated from isothermal, power-law, NFW and cored halos considered as prospective "universal" dSph halos (Section [5j. For each type 
of halo we fit only for the anisotropy and normalization. All isothermal, NFW and cored profiles above have normalization Vrnax ~ 10 — 20 km 
s^^ — see Table[3] All power-law profiles have normalization Afsoo ~ [0.5 — 1.5] X IO^A/q. 



by a and 7. Thus the parameter Ymax sets the normal- 
ization of the mass profile. 

The normalization can equivalently be set by specify- 
ing, rather than Vmax-, the enclosed mass at some par- 
ticular radius. For radius x, the enclosed mass M{x) 
specifies M(ro) according to 



M(ro) =M(x)- 



3— 7+Q . 



-1] 



rn / L a ' a 



3— 7+Q . 



(8) 

SOS demonstrate that for most dSphs the Jeans anal- 
ysis can tightly constrain Msoo- Here, in addition to 
A^sooj we shall consider the masses within two alterna- 
tive radii as free parameters with which to normalize the 
mass profile. Specifically, we consider the mass within 
the half-light radius, M{rhaif ), and the mass within the 
outermost data point of the empirical velocity dispersion 
profile, M{riast)- 



3.4. Markov-Chain Monte Carlo Method 

In order to evaluate a given halo model, we com- 
pare the projected velocity dispersion profile, ap{R), 
from Equation [3] to the empirical profile, (Tvb(i?), dis- 
played in Figure [T] For a given parameter set S = 
{— log(l — /3),logMx,logro,Q;,7}, where Mx is one of 
{Vmax, M{rhaif ), M^qq Or M{riast)}, we adopt uniform 
priors and consider the likelihood 



N 



exp 



1 iay.iRi) - apiR,)y 
'2 Var[(7y„(i?.)] 

(9) 

where Varfavb (^i)] is the square of the error associated 
with the empirical dispersion. 

Our mass models have five free parameters (four halo 
parameters plus one anisotropy parameter). In order 
to explore the large parameter space efficiently, we em- 
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ploy Markov-chain Monte Carlo (MCMC) techniques. 
That is, we use the standard Metropoli s-Hastings al- 
gorithm (|Metropolis et al.|[l953l : lHastingsill97Q) to gen- 
erate posterior distributions according to the following 
prescription: 1) from the current location in parame- 
ter space, Sn, draw a prospective new location, 5", from 
a Gaussian probability density centered on Sn', 2) eval- 
uate the ratio of likelihoods at Sn and 5"; and 3) if 
C{S')/({Sn) > 1, accept such that Sn+i = S', else ac- 
cept with probability C(5')/C(S'„), such that Sn+i = S' 
with probability (^{S')/C{Sn) and Sn+i = Sn with prob- 
ability 1-C(5')/C(5„). 

For this proce dure we use the ada ptive MCMC en- 
gine CosmoMC^ (jLewis fc Bridlell2002[ ). Although it was 
developed specifically for analysis of cosmic microwave 
background data, CosnioMC provides a generic sampler 
that continually updates the probability density accord- 
ing to the parameter covariances in order to optimize the 
acceptance rate. For each galaxy and parameterization 
we run four chains simultaneously, allowing each to pro- 
ceed until the variances of parameter values across the 
four chains become less than 1% of the mean of the vari- 
ances. Satisfaction of this convergence criterion typically 
requires ~ 10^ steps for our chains. We then identify the 
posterior distribution in parameter space with the last 
70% of all accepted points (we discard the first 30% of 
points, which correspond to the "burn-in" period.) 

3.5. Results 

Figure [2] indicates posterior distributions of model pa- 
rameters that we obtain from our MCMC analysis of 
Fornax, for each choice of parameter used to specify the 
normalization of the mass profile. We ran CosmoMC in- 
dependently for each normalization parameter, so while 
the distributions for other parameters should be similar 
for a given galaxy, they need not be identical. Contours 
in Figure [5] identify regions of 68% and 95% confidence. 

In general there exist degeneracies among the five pa- 
rameters such that most cannot be constrained uniquely. 
However, we find that most allowed mass profiles, regard- 
less of their shapes, intersect near the half-light radius. 
Then, as Figure [2] demonstrates, when we normalize the 
mass profile by either M{rhaif) or Afaoo, the degeneracies 
are such that the allowed region in parameter space sub- 
tends a relatively small mass range. Therefore while our 
analysis cannot place meaningful constraints on most pa- 
rameters, it does constrain Mivhaif) and M^qq. At much 
smaller and much larger radii (e.g., riast), the allowed 
mass profiles can diverge and the resulting constraints 
are weaker. 

In the interest of brevity we do not include the equiv- 
alent of Figure [5] for each of the other galaxies in our 
sample. Instead, for each of the eight dSphs with veloc- 
ity dispersion profiles in Figure [1] Figure [3] indicates the 
posterior distributions of M{rhaif) and anisotropy. Table 
[5] then lists the constraints on M{rhaif), Afaoo, M{riast) 
and V„iax- 

3.5.1. Mass 

We confirm the SOS result that the Jeans/MCMC anal- 
ysis places relatively tight constraints on M300, and that 

® available at [http : / / cosmologist . info / cosmomc | 



for the eight galaxies considered here, M300 ^ IO'^Mq. 
We obtain similarly tight constraints on M{rhaif), which 
is not surprising, since for these galaxies rhaif 300 pc. 
These results are not sensitive to the adopted form of the 
stellar density profile — we obtain similar values when we 
repeat the analysis using exponential and King profiles. 
Tablc[2]lists the masses and 68% confidence intervals ob- 
tained from the MCMC analysis. 

3.5.2. Vmax for Fornax 

Using simi la r me thods with smaller da t a sets , 
iStrigari et all (|2006[ ) and iPefiarrubia et all ()2008aD 
demonstrate a degeneracy between parameters Vmax and 
ro such that halos with arbitrarily large Vmax can, 
given the freedom to make also arbitrarily large, all 
have similar masses within the luminous region. How- 
ever, whereas most previously pub hshed velocity dis- 
persion pro files are quite fiat'' Ce.g.. iKlevna et al.ll2002l : 
IKoch et al.ll2007bllal : IWalker et al.lf200rbir we now detect 
gently decreasing dispersion in the outer regions of For- 
nax (Figure [1]). In our Jeans/MCMC analysis this be- 
havior constrains the scale radius of Fornax's dark mat- 
ter halo. Although the constraint is loose, ro ~ 0.49^0 45 
kpc, it is enough to help break the degeneracy with Vmax', 
in tests during which we artificially boost the velocity dis- 
persion in the outer bins such that Fornax's dispersion 
profile would be perfectly flat, our analysis fails to place 
upper limits on cither the scale radius or Vmax- Thus 
with the addition of the outer Fornax data, we now ob- 
tain the first model-independent constraint on Vmax for 
any dSph galaxy: Vmax — IS^s km s^^. For the remain- 
ing dSphs, the absence of any constraint on rg implies 
that we can place only lower limits on Vmax , typically of 
Vmax > 10 km s-1 (see Tabled]). 

One caveat regarding the apparent constraints on For- 
nax's Vmax and ro follows from the fact that we have 
assumed the stellar velocity anisotropy, f3, is constant. 
Our models therefore associate the onset of declining ve- 
locity dispersion with a transition in the shape of the 
density profile, which in our models helps to specify ro. 
In reality, falling dispersions may alternatively signal the 
onset of significant anisotropy. Therefore any constraints 
on ro, and in turn on Vmax, must be interpreted with 
some caution. However, the same is not true of the bulk 
mass constraints — e.g., SOS's models allow for radially 
variable anisotropy and demonstrate that constraints on 
M{rhaif) and M300 are robust over the full range of mod- 
els considered. 

3.5.3. No Constraint on Cores/Cusps 

For no dSph does our Jeans analysis place meaningful 
constraints on the halo sh ape parameters a and 7. Fo r 
the reasons explained by lEvans. An fc Walkeii (|2009f) . 
any such constraint on the inner density slope, 7, is 
strongly sensitive to the shape of the adopted stellar den- 
sity model*. Meaningful constraints on 7 require more 
sophisticated analyses that incorporate realistic models 

[Wilkinson et al.l 1I2OOI ') and IKlevna et all 1I2OOI ') report sharp 
declines at large radii in Ursa Minor and Sextans, but these fea- 
tures are not repro duced either in our data sets or in those of 
IMufioz et al.l ll2005t) . In any case such a sudden decline is difficult 
to reconcile with equilibrium models. 

* For a simple illus trat ion, consider the mass profile given by the 
first part of Equation 1101 which is derived from the Jeans equation 
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Fig. 2. — Two-dimensional posterior distributions of Fornax halo model parameters from Markov-Chain Monte Carlo (MCMC) analysis (section 
[31 ■ Contours identify regions of 68% and 95% confidence, respectively. Histograms indicate marginal distributions of each parameter. 



of the stellar distribution function fe.g.. [Wilkinson et alj 
l2002t Wilkinson et al. in prep.) and can thereby be 
evaluated on a star-by-star basis. 

3.6. A Simple and Robust Estimator for M(rhaif) 

We conclude that the available data tightly constrain 
dSph masses within r ^ r^aif- Furthermore, these 
constraints hold over a wide range of possible halo 
shapes and stellar velocity anisotropies. SOS reach the 
same conclusion regarding the robustness of M^qq, while 
IPefiarrubia et al. (2008a) show that M{rhaif) is well con- 
strained for an NFW halo regardless of the particular val- 
ues of Vmax and rp. Our analysis thus confirms the S08 
result a nd generalizes the conclusion of IPefiarrubia et all 
(|2008aD to include non-NFW halos. In short: we know 
M{rhaif) for the dwarf spheroidals. 

Because the MCMC method requires significant com- 

in the special case of /3 = and =const. Then the shape of 
M{r) is determined uniquely by the stellar density profile. 



putational effort, it is reasonable to wonder whether the 
measurement of such a bulk quantity as M{rhaif) might 
adequately be made using more conventional techniques. 
Therefore let us define a simple analytic model that re- 
lates Mijhaif) to observed quantities. Specifically, sup- 
pose the stellar component is distributed as a Plum- 
mer sphere with a velocity distribution that is isotropic 
(/? = 0) and has constant dispersion ay^{R) = v'^ = a^. 
Since these conditions are broadly consistent with all 
available dSph data^, this exercise amounts to consid- 
eration of a subset of the models not ruled out by the 
MCMC analysis discussed above. Then the Jeans equa- 

^ Velocity dispersion profiles are not yet available for the faintest 
dSphs, for which small samples provide measurements only of cen- 
tral velocity dispersions. Thus our assumption that the faintest 
dSphs have fiat velocity dispersion profiles remains untested. 
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tion gives 



M{r) = -■ 



Gv dr G[l + rVrLf] 



(10) 



where the last expression is specific to the assumption of 
a Plummer profile for the stellar density. Equation [TO] 
immediately provides the convenient estimate 



M{rhalf) = tJ-ThalfCr^ 



(11) 



where /i = 580 MQpc~^km~^s^, which differs only by a 
factor of 1.5 from the core-fitting formula com monly used 
to estima te dSph dynamical masses (e.g., IllingworthI 
[19.76; Ric hstone fc Tre mainc 1986; Mateo 199^^ 

Table [H lists estimates of M{rhaif) and M300 obtained 
from Equations [TOKTT] for the eight bright dSphs with 
velocity dispersion profiles, and Figure [3] indicates the 
estimates obtained from Equation [TT] for easy compar- 
ison to estimates obtained from the full Jeans/MCMC 



analysis. We find that the simple estimates from Equa- 
tion [11] generally stand in excellent agreement with con- 
straints we obtain from the full Jeans/MCMC analysis. 
Furthermore, as shown in Figure IH we reproduce the 
flat A/300-luminosity relation of S08 (c.f. their Figure 1) 
accurately merely by applying Equation [TO] to the data 
listed in Table [TJ We conclude that the mass estimates 
given by Equations 1 1 OlITT] are robust against a wide range 
of halo models when evaluated near the half-light radius. 
At radii much different from the half-light radius, the es- 
timates become model-dependent. In what follows, the 
values of Mirhaif) that we consider are those obtained 
from Equation [Til 

4. A NEW SCALING RELATION FOR DSPHS 

The left-hand panel of Figure [5] plots M{rhaif) (from 
Equation [TTjl against Vhaif for the 28 objects with pub- 
lished kinematic data (see Table [T]). There we find an ap- 
proximately power-law correlation between M(rhaif) and 
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TABLE 2 

Mass constraints from MCMC analysis and Eq. 10 



Galaxy 


Method 


A/(rhaif) 

r-1 n7 71 T 1 


Msoo 

ri r\7 71 /T 1 

[IO'A/q] 


riasf, M{riast) 

[kpc]; [1O'M0] 


n —11 
[km s 


Carina 
Carina 


MCMC 
Eq. 10 


— 0.1 

0.6 ±0.2 


0.6+»i 

— 0.2 

0.9 ±0.3 


0.87;3.7+?-i 


> 10 


Draco 
Draco 


MCMC 
Eq. 10 


— 0.3 

0.9 ±0.3 


— 0.6 

2.0 ±0.5 


0.92; 26.4^^7 4 


> 17 


Fornax 
Fornax 


MCMC 
Eq. 10 


4.3+0-? 

— 0.7 

5.3 ±0.9 


0.7+li 

— 0.2 

0.8 ±0.1 


1.7; 12.8+?-« 


18+5 

— 3 


Leo I 
Leo I 


MCMC 
Eq. 10 


1.0+0-^ 

— 0.4 
1.2 ± 0.4 


— 0.4 

1.8 ±0.5 


0.93; 8. 9+52 


> 15 


Leo II 
Leo II 


MCMC 
Eq. 10 


0.4 ±0.1 


-"-■■^-0.7 

1.2 ±0.3 


0.42;1.7Ii 2 


> 11 


Sculptor 
Sculptor 


MCMC 
Eq. 10 


1.3 ±0.4 


1 O+0.2 
^■''-0.2 

1.7 ±0.5 


1.1; lO.Otg Q 


> 15 


Sextans 
Sextans 


MCMC 
Eq. 10 


J^-"-0.4 

2.5 ± 0.9 


n 7+0.6 
-0.4, 
0.4 ±0.2 


1.0;2.0to ? 


> 9 


Ursa Minor 
Ursa Minor 


MCMC 
Eq. 10 


1 q+0.3 

^••^-0.5 

1.5 ± 0.4 


1 4+0.3 

1.7 ±0.4 


0.74; 4.4^2.0 


> 13 



we find log[crvo/(kms ^)] ^ 0.2 log[r,,a;//pc] + 0.3 (long- 
dashed line in Figure (5]). Using Equation [TT] to translate 
into the mass-radius plane, we obtain M{rhaif) / Mq ~ 
2300(r,ja///pc)^-^=^° '' (long-dashed line in Figurep. 

Constraints on M{rhaif) translate directly to give the 
mean density interior to the half-light radius, {p) = 
iM [r half) /{^T^T half)- The right-hand panel of Figure [5] 
plots mean density against half-light radius for the en- 
tire dSph sample. In terms of mean density inside Vhaif, 
the aforementioned relations imply (p)/(AfQpc~'^) ~ 
550(r/ia///pc)~^'^^° '*, and indeed, the data points in the 
right panel of Figure are well fit by this relation (long- 
dashed line). We note in particular that the dSph popu- 
lation adheres to this relation even as Vhaif reaches val- 
ues as small as tens of parsecs, with the caveat that in 
some cases the measured velocity dispersions, densities 
and masses are upper limits. If their velocity dispersions 
have indeed been resolved and reliably trace their masses, 
the smallest dSphs exhibit the largest galactic densities 
yet recorded, with {p) < SMqpc"^ ~ 200 GeV cm~^. 
These extreme values point to an early epoch of forma- 
tion, when the Universe itself would have had similar 
density, and may account for the survival of the ultra- 
faint systems against the destructive influence of tides 
(see below). 

4.1. The Effects of Tides 

The observable properties of the Local Group dSphs 
must be influenced to some degree by the tidal forces 
exerted by the M ilky Way and MB 1. Early work on 
this problem (e.g.. lOh et allll995l : iPTatek fc Prvodll995f) 
suggested that tides do not significantly alter the cen- 
tral velocity dispersion of a satellite as it passes through 
pericenter, lending credibility to estimates of dSph dy- 
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Fig. 4. — iV^soo? estimated using Eguation llOl versus to tal V-band lu- 
minosity. The flat relation reproduces the main result of lStrigari et all 
i200g|) (see their Figure 1). Filled circles represent galaxies for which 
observed kinematic tracers extend beyond r > 300 pc. Open squares 
represent ultra-faint satellites for which r^^lf < 100 pc, and for which 
there is no available evidence that a dark matter halo actually extends 
to a radius of 300 pc. 

fhaif- Of course, one expects power-law behavior merely 
from the form of Equation [Tl] — e.g., in the absence of 
a correlation between half-light radius and velocity dis- 
persion, one expects to fit the mass profile of a singu- 
lar isothermal sphere, M{r) cx r (dotted line in Figure 
[5]). Yet the slope in the empirical M{rhaif) — fhaif rela- 
tion is steeper than that corresponding to the isothermal 
sphere. This result follows entirely from the fact that, 
given the data now available for ultra-faint dSphs, there 
exists a clear correlation, shown in Figure [51 between 
dSph velocity dispersion and half-light radius. Fitting a 
power law in the plane of directly observed quantities, 
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Fig. 5. — Left: M{rhatf) against half-light radius for all measured dSphs and dSph candidates (see Tablc0. Over-plotted are the best-fitting 
mass profiles, M(r), for isothermal, power-law, NFW {Vmax — 15 km s~^, ro — 795 pc) and cored (Vmax — 13 km s~^, ro — 150 pc) halo profiles. 
Right: Mean density within the half-light radius, versus rf^^lf- Over-plotted are the same halo models as in the left panel. In both panels, arrows 
indicate the m agnitude and direction that individual galaxies would be displaced due to the tidal stripping of 99% of the original stellar mass 
dPenarrubia et al.u2008bO . Text markers identify the most extreme outliers and Sagittarius. 

of analytic formulae to describe the tidal evolution of 
the parameters relevant to the present work. In par- 
ticular, they find that for a dSph consisting of a stel- 
lar component described by a King profile embedded 
within an NFW halo, all parameters evolve according 
to the fraction of mass lost within the core radius (here 
taken to be approximately the half-light radius). For 
tidal stripping that results in the loss of as much as 99% 
of the original luminosity, the evolution of parameters 
is described by L/Lq « 2^(T/aQ)^-'^[l + crVcrg]-? ^^^^ 
Thaif/rhaifo « 23/2(ct/(To)1-3[1 + a^a^]'^/^, where sub- 
scripts of naught indicate initial values (P08). 

Arrows in Figures [5] and [5] indicate the directions and 
magnitudes of the tidal tracks followed by a satellite as it 
loses 99% of its stellar mass under the P08 relations. In 
the a-rhaif plane of obvervables (Figure [6]), at all points 
along the tidal track the displacement vector points pri- 
marily downward such that, as the satellite loses mass, 
its velocity dispersion decreases more quickly than its 
size. As a result, tides tend to carry individual objects 
off the empirical relation rather than move them along it. 
While the tidal track in the M{rhaif)-rhaif plane (Fig- 
ure [5l left) is more parallel to the empirical relation, the 
tidal track in the {p)-rhai / plane (right panel of Figure 
[5|) provides strong evidence against a tidal origin for the 
empirical relation. There we find that tides tend to de- 
crease both the half-light radius and the mean density 
within the half-light radius. Thus tides do not generate 
the circumstance in which the smallest objects also have 
the largest mean densities interior to rhai / ■ These consid- 
erations lead us to conclude that tides did not introduce 
the correlations exhibited in Figures [H] and [SI This is 
not to say that tides are irrelevant — we now discuss two 
ways in which tides may have altered the character of 
the correlations we now observe. 
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Fig. 6. — Global velocity dispersion against half-light radius for all 
measured dSphs and dSph candidates (see Table [3. Over-plotted are 
the scaling relations that correspond to the best-fitting mass profiles 
shown in Figure[5] For the NFW and cored halos, these arc translated 
into the (J-Vhaif plane via Equations 1 12|13l As in Figure [s] the arrow 
indicates the magnitude and direction of dis placement due to the tida l 
stripping of 99% of the original stellar mass llPeiiarrubia et al.|[2008bl) . 

namical masses. More recent N-body simulations (e.g., 
Read et al. 2006; Muiioz et al. 200^; iPenarrubia et all 
2008b; .Klimentowski et al.i ,2008 ) monitor the evolution 



of both stellar and dark matter components over sev- 
eral peric entric passages and a wide range of orbits. 
Feiiarrubia et al.l (|2008bl "P08" hereafter) provide a set 
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First, we can expect at least some of the scatter in the 
correlations to be due to the varying degrees to which 
dSphs on different orbits have been tidally sculpted. In 
this regard, the dSphs Cetus and Tucana merit specific 
attention because they are the only objects in our sam- 
ple that have no clear association with either the Milky 
Way or M31. Cetus lies at a distance of D ^ 755 kpc 
(|McConnachie fc IrwinI f2006f) and Tucana i s even more 
remote, at D ~ 880 kpc (jSaviane et al.lll996l ). It is there- 
fore reasonable to assume that neither object has lost ap- 
preciable amounts of mass to tidal stripping. Then it is 
perhaps not a coincidence that Cetus and Tucana are the 
two most extreme outliers with respect to the empirical 
relations in Figures [5] and [6] Both objects have large ve- 
locity dispersions (and hence large M{rhaif)) compared 
to other objects of similar size. We speculate that tides 
may have influenced all other dSphs to a significantly 
higher degree than Cetus and Tucana, moving the popu- 
lation systematically along POS's tidal tracks while leav- 
ing Cetus and Tucana in relatively pristine condition. 
Notice that even strong tidal stripping, if it acts to a 
similar extent on most dSphs, can preserve the slope of 
an initial correlation despite shifting the locus of the af- 
fected population. 

Second, tides might alter the slopes of the correlations 
in Figures [5] and [S] if systems of a particular size happen 
to be the most prone to tidal stripping. According to the 
tidal tracks, tides act to move objects primarily down- 
ward, towards lower velocity dispersion in the a-rhaif 
plane of Figure [HI But we must consider that since the 
objects with smallest half-light radii tend also to be found 
closest to the Milky Way, tides may have displaced them 
systematically farther along the tidal track than larger, 
more distant systems. This may have increased the slope 
of the relation in Figure[51 Under this scenario, the dSph 
population may have formed according to a flatter rela- 
tion better described by an isothermal sphere {M oc r), 
but then evolved tidally such that the relation we observe 
today is the steeper power law M oc r^ '^. However, we 
re-emphasize that while tides may have increased both 
the scatter and slope of the empirical correlations, move- 
ment of individual systems along tidal tracks would not 
have introduced correlations where none existed initially. 

4.2. Sagittarius 

Sagittarius (Sgr) is the only one of the Local 
Group dSphs that is unambi guous l y in the thr oes of 
tidal disruption (llbata et all Il995l : "iMateo et all 119961 : 
iMaiewski et al.ll2003[ ). The largest of the objects in our 
sample, the main body of Sgr has a half-light radius 
that is slightly larger than the scale radius of the best- 
fitting universal NFW halo. This is compatible with 
the fact that Sagittarius continues to lose stars, as its 
streams of debris wrap more than once around the sky 
(jMajewski et al. 2003). Perhaps strangely, Sgr is other- 
wise inconspicuous in our scaling relations, falling neatly 
onto the best-fitting power-law profile. According to 
P08's tidal tracks, the progenitor of Sgr had larger veloc- 
ity dispersion than Sgr exhibits today, such that Sgr may 
originally have been a more extreme outlier than Cetus 
and Tucana are now. 

4.3. On the Ultra-Faint Satellites 



Since they extend the range of dSph sizes downward 
by an order of magnitude, the ultra-faint dSphs provide 
much of the leverage in discerning the scaling relations 
underlying Figures [MSI Therefore we must note some 
caveats due to the fact that the least luminous objects 
are necessarily the least well-studied. First, their veloc- 
ity dispersions are the most uncertain, a result of small 
sample sizes (e.g., ju st five stars in the case of Leo V 
(| Walker et al.ll2009al) ') and the fact that their small dis- 
persions are near the resolutio n limits o f the spectro- 
graphs used to measure them (S imon fc G eha 2007). In 
some cases, the measured dispersions (and corresponding 
masses and densities) represent just upper limits, so the 
data for objects like Leo V and Bootes II remain con- 
sistent with these objects having tiny dispersions that 
would render them more similar to star clusters than to 
bona fide dSph galaxies. 

Second, many (e.g., BooII, Coma, Segue 1, Segue 2, 
UMall, Willman 1) of the ultra-faint satellites are found 
within < 50 kpc of the Sun, where they are most vul- 
nerable to Milky Way tides (see also section 14. 1[) . While 
extreme central densities (< 5M0pc~'^) should be suffi- 
cient to protect the smallest systems against tidal disrup- 
tion, many of the ultra-faint sate llites appear elongate d 
(e.g., Hercules has axis ratio 3:1 (jColeman et al.l [20071 )) 
or otherwise exhibit irregular morphologies (e.g., U Mall 
has an apparent double core ([Zucker et al.ll2006aD V On 
one hand, one does not expect unbound tidal debris to 
inflate the veloci ty dispersions of the smallest dSphs — 
iPenarrubia et al] ((2008) use simulations to show that 
unbound debris escapes from a disrupting dSph on a 
timescale similar to the crossing time, which for the 
smallest satellites is just ^ rhaif ^ 10 Myr. On the 
other hand, if distorted morphologies result from strong 
tidal processes, then it is difficult to understand how 
these systems could be embedded in intact dark matter 
halos similar to those inferred for more regular dSphs. 
It may turn out, then, that some of the recently discov- 
ered satellites are unbound star clusters in the process of 
dissolving. In that case, it would be puzzling that these 
systems manage to fall on or near the empirical relations 
defined by bona fide dSphs. 

5. A UNIVERSAL MASS PROFILE? 

We have shown that there exists in the population 
of known dSphs an empirical correlation between half- 
light radius and velocity dispersion (Figure [S]) that, un- 
der equilibrium models, implies a correlation of the form 
M{riiaif ) oc rj^'^j^'^'* (Figure [5]). We emphasize that this 
implication is insensitive to the assumption of isotropy 
inherent in Equations [TUKTTl — our Jeans /MCMC analy- 
sis shows that our estimates of M{rhaif) are robust over 
a wide range of (constant) anisotropy values. Taken at 
face value, these results suggest that the smallest dSphs 
may simply be embedded more deeply inside dark mat- 
ter halos that are otherwise similar to those inhabited 
by the larger (and more luminous) dSphs. Then let us 
examine the extent to which dSph dark matter halos 
might be "similar." In fact let us carry this notion to 
its limit and consider the hypothesis that all dSphs are 
embedded within identical dark matter halos. In that 
case our robust constraints on M(rhaif) become data 
points that sample a single, "universal" mass profile at 
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discrete radii, thereby allowing us to measure the pro- 
file's shape directly and free from the concerns regarding 
mass/anisotropy degeneracy that plague analyses of indi- 
vidual systems. We have already shown that the data are 
broadly consistent with a power-law profile of the form 



M{r) 



(X r 



In what follows we also consider cusped 



and cored profiles, and then we examine the ability of a 
universal profile to fit the individual velocity dispersion 
profiles of the brightest dSphs. Finally, we quantify the 
amount of empirical scatter with respect to each candi- 
date for a universal dSph profile. 

5.1. Cusps versus Cores 

We now investigate whether the correlations in Figures 
[5] and [S] can be described in terms of universal versions of 
the halo models considered in Section [31 For simplicity 
we restrict our consideration to two specific halos of in- 
terest: the cuspy NFW halo with a — j = 1 and a cored 
halo with a = 1, 7 = 0. In order to compare these halo 
models directly with observable quantities, we equate 
M[rhaif) from Equation [7] with the model-independent 
estimate given by Equation [TT] For the NFW halo, this 
gives the scaling relation 



5 ln[l + ^]-^ 



(12) 



where 77 ~ 2.16. For the cored halo, the scaling relation 
becomes 



rhaifcr 



2r]roV- 



7nax 



5{HI + V] + T^-2(TTW-I^ 
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now with r] ~ 4.42 for a = 1, 7 = 0. 

Let us consider whether the correlations in Figures [5] 
and [6] can be fit with a single halo of either the cusped or 
cored variety. Contours in the Vmax-ro plane in Figure [7] 
indicate constraints we obtain by fitting, via Equations 
fT2landfT3l NFW and cored profiles to the rhaif-ava data. 
Accounting for measurement errors in both dimensions, 
the best- fitting NFW halo has Vmax ~ 15 km s~^, ro ~ 
795 pc, while the best-fitting cored halo has Vmax ^ 13 
km s~^, ro ~ 150 pc. The best-fitting NFW and cored 
profiles are plotted over the empirical relations in Figures 
El and El 

For both NFW and cored halos the scale radius is 
constrained by the slope in the empirical cr-rhaif rela- 
tion (Figure [6]). Cored halos with scale radii larger than 
ro > 200 pc generate correlations with slopes in the a- 
fhaif plane that are systematically steeper than that of 
the empirical relation. We have tried other forms of the 
cored halo (e.g., a — 0.5, a = 2) and find that in all 
cases, the scale radius must be < 200 pc. Another way 
of understanding this constraint is that it is set by the 
fact that the empirical relation in Figure El (right panel) 
does not show signs of turning over as r — > 0. Thus the 
universal dSph halo, if it is cored, must have small scale 
radius. 

Figure 13 also demonstrates that for a universal 
NFW halo, the allowed region of parameter space in- 
tersects the mass-concentration relation that is de- 
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Fig. 7. — Constraints on halo parameters Vmax and rp if we fit a 
single NFW (blue) or cored (red) halo to the vs. Vhaif relation for 
all dSphs. Contours indicate 68% and 95% confidence regions, from a 
fit. For comparison, solid black contours show model-independent 
constraints obtained for the Fornax halo alone (Section |3.4| and Fig- 
ure [2}. Overplotted are the mass-concentration relations for cold dark 
matter halos at z — 0,1,2, f rom cosmolo^ical N-body simulatio ns 
INavarro. Frenk fc Whitgil996l : lEke et al.ll200ll : IBullock et"alll2001l) . 

rived i ndependently from cosmological N-body simu- 
lations (INavarro. Frenk fc Whitd 119961: l Eke et all 120011 : 
IBullock et al.l 120011 : iPefiarrubia e t al. 2008aD Thus the 
single NFW halo that best fits all the data is broadly 
consistent with the ACDM framework. 

From both the cusped and cored halo fits there is a 
clear prediction that deviates from that of the power-law 
fit: collectively, dSph velocity dispersions increase with 
half-light radius until reaching a maximum of cr '--^ 10 km 
s^^ for r^aif ~ 1 kpc (Figure E]). Larger systems falling 
on the same relation will have smaller velocity disper- 
sions, cr < 10 km s^^, as the half-light radius becomes 
similar to the scale radius of the dark matter halo. We 
will soon be able to test this prediction as kinematic data 
become available for more dSph s atellites of M31, many 



of which have 



rhalf 



> 



1 kpc (jMcConnachie fc Irwin! 
|2006( ). The universality of the best-fitting cusped or 
cored halos would imply that dSphs represent a class 
of spheroid that is at once distinct from globular clus- 
ters {rhalf 1 — lOpc and cr ~ 5 — 10 km s~^) and 
larger dwarf elliptical galaxies {rhalf ~ 1 kpc, cvb ~ 50 
km s"""^) — both types of object have velocity dispersions 
too large for the dSph relation for cored or cusped halos. 
In contrast, a power-law profile would imply that veloc- 
ity dispersion increases with half-light radius indefinitely 
(dashed line in Figure E]), admitting the possibility that 
dSph mass profiles are similar to those of larger ellipticals 
and perhaps spiral galaxies (Paulo Salucci, private com- 
munication; Stacy McGaugh, private communication). 

Irrespective of whether the halo is cored or cusped, the 
hypothesis of universality leads to another clear predic- 
tion. If the right-hand side of Equation [1] is universal, 
then the left-hand side must be as well. So, there is a 
universal functional relationship between the dSph stel- 
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lar density, velocity dispersion and anisotropy. This is 
simple to work out in the case of isotropy and constant 
velocity dispersion. Two dSphs with velocity dispersions 
fJi and CT2 must have luminosity profiles vi and 1^2 related 

by 

= (14) 

In other words, the two profiles are related via a power- 
law transformation. Even in the more complicated case 
when the anisotropy /3 is non-zero, the luminosity pro- 
files of different dSphs are still related to each other, 
though the transformation is now more complicated than 
a power-law. 

5.2. Velocity Dispersion Profiles 

A universal dSph dark matter halo must account not 
only for the bulk kinematic properties of the dSph pop- 
ulation, but also for the velocity dispersion profiles, 
where available, of individual dSphs. Therefore let us 
test whether the four candidate profiles considered in 
this work — a singular isothermal sphere with M oc r, 
a power-law profile with M (x r^-*, an NFW cusp with 
scale radius rg = 795 pc, and a core (a = 1; 7 = 0) with 
ro — 150 pc — are consistent with the velocity dispersion 
profiles in Figure [T] For each profile we fit to the empir- 
ical velocity dispersion profiles using the Jeans equation 
(Equation[3|) , allowing now for just two free parameters — 
anisotropy and the normalization, Vmax- For the power- 
law profile the circular velocity curve is unbounded as 
r ^ 00, so we normalize the power-law profile with the 
value of M300. In order to consider only the most real- 
istic cases, we restrict the anisotropy to values between 
-0.6 <I3< -hO.3. 

For each of the four universal profiles considered, Table 
[3] lists the values of (3 and Vmax that produce the best fits 
to the empirical velocity dispersion profiles of the bright 
dSphs. We find that for each halo profile, the data for all 
eight dSphs are well fit if we allow the normalization to 
vary over a factor of < 2. For the power-law profile, we 
can fit all the empirical velocity dispersion profiles if we 
let M300 take values between [0.5 — 1.5] x l(f Mq — note 
that this is the same range for M300 found by SOS, but 
unlike in that study, here M300 serves as the normaliza- 
tion for mass profiles of the same shape. The isothermal, 
cored and NFW profiles require Vmax in the range 10 — 20 
km s^^. The best fits for each profile are plotted against 
the empirical velocity dispersion profiles in Figure [U and 
parameters of those fits are listed in Table [3l 

5.3. Scatter 

We now evaluate the feasibility of a universal dSph 
mass profile in terms of the amount of scatter with 
respect to the best candidates for that profile. Fig- 
ure [5] plots residuals A(logM) = \og[M{rhaif)lMQ\ - 
\og[M{rhaif)universailMQ\ (where the first term rep- 
resents the data and the second term represents 
the universal profile) with respect to the best-fitting 
isothermal, power-law, NFW and cored profiles con- 
sidered in this work. For comparison the bottom 
panel plots the residuals in the common-Msoo re- 
lation, log[M3oo/M0] — 7, where M300 is estimated 
from Equation [TOl Text indicates for each profile 

the statistic = E^=l(A(logM))V<5^ where 



5 = <JM(rhaif )[^^{'''haif ) In(lO)] represents the measure- 
ment error in log-space. The power-law [Mir) cx r^ **) 
provides the best fit (x^ = 2.8), followed by the NFW 
profile (x^ = 3.9), the cored profile (x^ = 5.3), and the 
isothermal profile (x^ = 12.6). While these values of x^ 
indicate scatter larger than what is expected from the 
quoted errors if the profile is truly universal, we note 
that that the scatter of Afsoo values about a common 
mass of lO^Af0 corresponds to x^ = 10.0. In addition 
to x^ J Figure [5] also indicates the root-mean-square of 
the residuals. There we also find that the scatter about 
a universal mass profile (rms — 0.33 for the power-law 
profile) is similar to the scatter about a common value 
of A%o {rms = 0.41). 

Thus in terms of empirical scatter, the hypothesis that 
dSphs follow a universal mass profile receives as much 
support from the data as the claim that all dSphs have 
a common M^qq. Note that the former is the stronger 
claim, since it implies the latter (not vice- versa), and 
both statements involve an extrapolation of sorts. SOS's 
argument for a common M300 rests on extrapolation, 
since for the smallest objects (ruai / ~ 30 pc) it requires 
evaluation of mass profiles at radii that are larger by an 
order of magnitude than the region containing kinematic 
data. In contrast the claim for a universal mass profile 
emerges from the pattern of masses robustly measured 
at the half-light radius of each dSph. Yet it is important 
to acknowledge that we know the enclosed mass only 
at this one particular radius in each dSph, so we can- 
not definitively rule out the situation in which any given 
dSph has a mass profile that diverges wildly from the 
hypothesized universal profile at radii far from its own 
half-light radius. However, if this circumstance were true 
for more than just a few outliers, then it would require a 
remarkable coincidence for the masses at dSph half-light 
radii to generate the appearance of a universal profile 
over two orders of magnitude in rhai f ■ We conclude that 
if one is impressed by the commonality of Af3oo among 
dSphs, one is obliged to be impressed by the apparent 
commonality of dSph masses at all radii. 

6. DISCUSSION AND SUMMARY 

In the first part of our analysis (Section [3]), we have 
applied the Jeans equation to the velocity dispersion 
profiles of eight bright dSphs. The quality of the data 
set now available for Fornax allows us to place the first 
model-independent constraint on the maximum circular 
velocity of a dwarf spheroidal. While for all other galax- 
ies we obtain only lower limits of Vmax ^ 10 km s^^, for 
Fornax our Jeans/MCMC anlysis yields Vmax ~ 18^3 
km s^^. This constraint enables straightforward com- 
parison to subhalos produced in cosmological N-body 
simulations, for which Vmax is readily quantifiable, and 
motivates the pursuit of larger velocity samples in the 
remaining dSphs. 

A key result from our Jeans/MCMC analysis is that 
for all well-observed dSph galaxies with measured ve- 
locity dispersion profiles, most allowed mass profiles in- 
tersect near the half-light radius, regardless of the as- 
sumed aniso tropy. This result is c onsist ent with previous 
finding s by iStrigari et al.l ()2008[ ) and iPeiiarrubia et all 
(|2008al) . and has two important consequences. First, 
while the Jeans analysis does not place meaningful con- 
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TABLE 3 

Anisotropy and Normalization of Universal Profile, for dSphs in Figure 1 



M ocr M PC r^-'^ NFW (ro = 795pc) Core (ro = 150pc) 

Galaxy /3; Vmax P; M30O P; Vmax P\ Vmax 

[Mq] [Mq] [km s-l] [km s-i] 

Carina -0.6+0-5; loll -0.6+0-9; oj+0.2 ^ -OAt^of, ISlJ -0.5±0;J; lit? 

Draco -0.5l°;5; ISt^ -0.6+0-9; i ^ ^^7 _o.6+0-8; 20+1 -0.6+0-^ 17tl 

Fornax -0.6+0-2; ^f^+i _o.2lO-3. 1.2+O.1 x 10^ -0.310-3. ig+i -O.SlO-l; 16^2 

Leo I -0.6+0-7; 1612 _o.6+0-9; l.e^O-* x 10^ -0.210-6. aol^ -0.6+0-9; 18+2 

Leo II -o.4lO;^; 12I2 +o.3_o.9; i-2lo;5 x lo'^ +o.3_o.9; I8I3 +o.3_o.9; 15I3 

Sculptor -0.6+0-2; 14+1 _o.5lO-3; i.il0.4 x lo^ -0.5lO-j; 19+\ -0.6+0-4; ig+l 

Sextans -O.SlO;?; 14t\ -0.510-J; 0.5l0;2 x 10^ -0.410-J; I1I2 -0.5lO;«; lol^ 

UMi -0.6+0-5; 14+2 _o.6+o.7. 1.3+0.4 X io7 -0.6+0-8; 19+2 _o.6+0-7; 1512 



straints on most halo parameters for most galaxies, for 
all eight "classical" dSphs we obtain relatively tight 
constraints on the enclosed mass at the half-light ra- 
dius. Second, the insensitivity of M(rhaif) to anisotropy 
means that one need not employ the full Jeans/MCMC 
analysis to obtain reliable estimates of M(rhaif )- In Sec- 
tion 13.61 we derive Equations [10] and [Tl] from the sim- 
plifying assumptions that dSph stellar components have 
isotropic velocity distributions -with constant dispersion. 
These mass estimators have the advantage that they re- 
quire as input only measurements of the velocity dis- 
persion and half-light radius, and so can be applied to 
the faint dSphs for which velocity dispersion profiles are 
unavailable. Using Equation (TU] to estimate M300, we 
immediately reproduce the flat Afaoo-himinosity relation 
(FigureU]) originally obtained by the Jeans/MCMC anal- 
ysis of SOS. Using Equation [HI we obtain estimates of 
M{rhaif) that stand in excellent agreement with the re- 
sults of our own Jeans/MCMC analysis (Figure [3] and 
Table O for the brightest dSphs. 

In the second part of this work, we have applied the 
simple estimator of M{rhaif) given by EquationfTTlto the 
entire dSph population in lieu of the full Jeans/MCMC 
analysis. We find that the 28 dSphs with published kine- 
matic data exhibit a correlation of the form M(rhaif) oc 
r^'^j^"'^ (Figure [51 left panel). This slope is steeper than 
that of the singular isothermal sphere because, as the 
inclusion of the ultra-faint dSphs makes clear, veloc- 
ity dispersion increases with half-light radius (Figure [6|). 
Constraints on M(rhaif) imply constraints on the mean 
density interior to the half-light radius, and the small- 
est dSphs now allow us to measure galactic densities at 
the scale of tens of parsecs. Toward these small radii, 
the mean densities we measure continue to rise accord- 
ing to the power law that best fits the empirical rela- 
tion, (p) (X r^^-^^OA (pigure [sj right-hand panel). For 
the smallest systems, this behavior may imply extremely 
large densities of order (p) ~ 5Mopc~ ^ 200Gev cm~^. 
However, we emphasize that for some of the smallest sys- 



tems, the available data do not place firm lower limits on 
velocity dispersion, so these densities should be viewed 
as upper limits. 

We have argued in Section 14.11 that tidal forces have 
not introduced the scaling relations underlying Figures[5]- 
[6l According to^the tidal tra cks p rovid ed by the N-body 
simulations of jPeharrubia et al.l (j2008bh . the long-term 
effect of tidal stripping is to reduce a dSph's velocity dis- 
persion faster than its size. This moves a dSph off the em- 
pirical fJ-rhaif relation rather than along it. Thus tides 
act primarily to increase the scatter in the scaling rela- 
tions we observe today. We conclude that the correlation 
between velocity dispersion and r^aif, and the resulting 
relation between M{rhaif) and ruaif, are inherent to the 
dSph formation mechanism. However, tides may have 
altered the slope of the empirical correlations we observe 
today. Since the dSphs closest to the Milky Way tend 
also to be the smallest and least luminous, the stronger 
tides they have encountered may have shifted their ve- 
locity dispersions systematically downward, thereby in- 
creasing the slope in the cr-rhai / relation. 

Finally, in Section [5] we have considered, in light of the 
correlation between M[rhaif) and rhaif, the possibility 
that all dSphs follow a universal mass profile. We have 
shown that in addition to the simple power law M(r) oc 
r^-^, the empirical relationship between mass and half- 
light radius (or equivalently, between velocity dispersion 
and half-light radius) can be fit by an NFW {Vmax = 15 
km s~^, ro = 795 pc) or cored {Vmax = 13 km s~^, 
ro = 150 pc) halo profile. The allowed parameters of 
the universal NFW halo (Figure [7]) overlap the mass- 
concentration relation derived from cosmological N-body 
simulations, supporting the internal consistency of the 
ACDM framework. While we obtain the best fits to the 
velocity dispersion profiles of the brightest dSphs if we 
allow the normalizations of the profiles to vary over a 
factor of < 2 (Section 15. 2p . we find that the empirical 
scatter of the M{rhaif) — rhaif data about a universal 
profile is similar to the scatter about a common value 
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Fig. 8. — Top four panels: Residuals with respect to universal mass 
profiles considered in this work. Bottom: Residuals with respect to the 
common mass scale -A/300 — 10^A/q. Empirical masses are estimated 
from Equations IIOIIII In the bottom panel, open squares represent 
ultra-faint satellites for which ri^^lf < 100 pc, and for which there is 
no available evidence that a dark matter halo actually extends to a 
radius of 300 pc. 



of A/300 • Thus the stronger claim that dSphs follow a 
universal mass profile receives at least as much support 
from the data as the independent claim that all dSphs 
have a common Maoo- 

It has long been suggested that dSphs exhibit a com- 
mon mass scale. Fifteen years ago, this conclusion fol- 
lowed simply from the virial theorem and the fact that 
the bright dSphs known at the time {Ly ^ 10^Lv,(7,) all 
have similar sizes and velocity dispersions (jMateo et al.l 
119931) . The recently-discovered ultra-faint dSphs, which 
are systematically smaller and kinematically colder than 
the brighter dSphs, undermine a com mon dSph mass 
scale as derived from the virial theorem. iSimon fc Gehal 
(j2007l ) show that the ultra-faint dSphs deviate from the 
luminosity-M/L relation followed by the brighter dSphs 
(Mateo 199^), and the dynamical masses of the small- 
est d Sphs ar e just - IO'^Mq (e.g., IGeha et all I2009bl : 
.Walker et all SOS are able to recover a com- 

mon dSph mass scale of M300 ~ lO'^Af© that includes 
all ultra-faint as well as bright dSphs, but doing so re- 
quires them to evaluate all mass profiles at a fixed radius 
that happens to lie well beyond the stellar component 
in the smallest dSphs. In contrast, we have presented 
the case for a common dSph mass profile that emerges 
from masses that relate directly to empirical dynami- 
cal properties — tracer density and velocity dispersion — 
in the regions where they are directly constrained by 
data, requiring no extrapolation to radii beyond the op- 
tical extent of the galaxy. 
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